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I i' ABSTRACT 

^ : 

Qh' We investigate the effects of clianges in the cosmological parameters between the 

Q I WMAP 1st, 3rd, and 5th year results on the structure of dark matter haloes. We 

5^ ■ use a set of simulations that cover 5 decades in halo mass ranging from the scales 

' of dwarf galaxies (Vc ~ 30 km/s) to clusters of galaxies (Vc ~ 1000 km/s). We find 

, that the concentration mass relation is a power law in all three cosmologies. However 

the slope is shallower and the zero point is lower moving from WMAPl to WMAP5 
^ ; to WMAP3. For haloes of mass logM2oo/[/i"^ Mq] = 10, 12, and 14 the differences 

, in the concentration parameter between WMAPl and WMAP3 are a factor of 1.55, 

\^ • 1.41, and 1.29, respectively. As we show, this brings the central densities of dark 

I matter haloes in good agreement with the central densities of dwarf and low surface 

Q> i brightness galaxies inferred from their rotation curves, for both the WMAP3 and 

7—i ' WMAP5 cosmologies. We also show that none of the existing toy models for the 

ly-^ , concentration-mass relation can reproduce our simulation results over the entire range 

■ of masses probed. In particular, the model of Bullock et al. (2001a; hereafter BOl) 

I fails at the higher mass end {M > lO^'^h^^ M©), while the NEW model of Navarro, 

O ■ Frenk & White (1997) fails dramatically at the low mass end (A/ < lO^^h''^ Mq). We 

^ ' present a new model, based on a simple modification of that of BOl, which reproduces 

. ^ , the concentration-mass relations in our simulations over the entire range of masses 

^ ■ probed (lO^^'/i-i Mg < M < lO^^h-iMo). Haloes in the WMAP3 cosmology (at a 

?— ( ' fixed mass) are more flatted compared to the WMAPl cosmology, with a medium to 

long axis ration reduced by « 10%. Finally, we show that the distribution of halo spin 
parameters is the same for all three cosmologies. 

Key words: galaxies: haloes - cosmology:theory, dark matter, gravitation - methods: 
numerical, N-body simulation 



1 INTRODUCTION 

In the paradigm of hierarchical structure formation, dark 
matter (DM) haloes provide the potential well in which 
galaxies subsequently form. As a consequence the structural 
parameters of disk galaxies (size and rotation velocity) are 
tightly coupled with those of their hosting DM halo, such as 
concentration and spin (e.g. Mo, Mao & White 1998; Dutton 
et al. 2007). 

It has been shown by several studies that the struc- 
tural properties of dark matter haloes are dependent on 
halo mass: for example higher mass halos are less concen- 
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trated (Navarro, Frank & White 1997, hereafter NFW; Bul- 
lock et al. 2001a, hereafter BOl; Eke Navarro & Steinmetz 
2001; Kuhlen et al. 2005; Neto et al. 2007; Maccio et al. 2007, 
hereafter M07), and are more prolate (Jing & Suto 2002; All- 
good et al. 2006; M07) on average. In the case of the spin 
parameter, however, there seems to be no mass dependence. 

In M07 we used a set of numerical simulations of struc- 
ture formation in a ACDM cosmology, with cosmological 
parameters that were motivated by the first year results of 
the Wilkinson Microwave Anisotropy Probe (WMAP) mis- 
sion (Spergel et al. 2003; hereafter WMAPl), to study how 
concentrations, shapes and spin parameters of dark matter 
haloes scale with halo mass. In this paper we extent the 
M07 study by investigating how these scaling relations de- 
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Table 1. Cosmological Parameters 



Table 2. N-body Simulation Parameters 



Name 






h 




n 




WMAPl 


1.0 


0.268 


0.71 


0.90 


1.00 


0.044 


WMAP3 


1.0 


0.238 


0.73 


0.75 


0.95 


0.042 


WMAP5 


1.0 


0.258 


0.72 


0.796 


0.963 


0.0438 



pend on the adopted cosmology. In particular, we present a 
large suite of numerical simulations for ACDM cosmologies 
whose parameters are motivated by the three and five year 
data releases of the WMAP mission (Spergel et al. 2007; Ko- 
matsu et al. 2008). In what follows we refer to these cosmolo- 
gies as the WMAP3 and WMAP5 cosmologies, respectively. 
The parameters of the WMAPl, WMAP3 and WMAP5 
cosmologies are listed in Table [T] Note that the WMAP3 
and WMAP5 cosmologies both have a lower matter density, 
Q.m, and a lower power spectrum normalization, erg (defined 
as the rms of the matter field, linearly extrapolated to the 
present, on a scale of 8/i~^Mpc) than the WMAPl cosmol- 
ogy. This implies that dark matter haloes of a given mass 
assemble later (e.g., van den Bosch 2001). Since the con- 
centration of a dark matter halo is related to its assembly 
redshift (e.g., Wechsler et al. 2002; Zhao et al. 2003a,b; Li 
et al. 2007), dark matter haloes are expected to be less con- 
centrated in the WMAP3 and WMAP5 cosmologies. This 
may help to reconcile the conflict between the concentration 
parameters inferred from the rotation curves of dwarf and 
low surface brightness (LSB) galaxies (de Blok, McGaugh & 
Rubin 2001; Swaters et al. 2003) and those predicted from 
the ACDM scenario. In addition, a lower amplitude of the 
power spectrum results in haloes being more prolate (All- 
good et al. 2006), which may have observational implications 
for the rotation curves of dark matter dominated galaxies. 

As in M07 we use a large suite of N-body simula- 
tions for the WMAPl, WMAP3 and WMAP5 cosmolo- 
gies with different box sizes to cover the entire halo mass 
range from 10^° M0 (haloes that host dwarf galaxies) 
to lQ^^h~^ M0 (massive clusters). We use these simulations 
to investigate the concentrations, spin parameters of shapes 
of dark matter haloes. We also present a critical compari- 
son of the simulation results with two toy models for the 
concentration-mass relation. In particular, we show that 
the BOl model is only able to reproduce the simulation re- 
sults in the galaxy mass range, but underpredicts the con- 
centrations of massive haloes. On the contrary, the NFW 
model yields a reasonable fit at the massive end, but dra- 
matically underestimates the concentrations for haloes with 
M < lO^^/i"^M0. We present a new model, based on a 
simple modification of the BOl model, that that is able to 
reproduce the concentration-mass relation over the entire 
mass range probed by our numerical simulations. 

This paper is organized as follows: in Section[2]the sim- 
ulations and the determination of the halo parameters are 
presented. In Sections |31 H] and [S] we discuss the results for 
halo concentrations, spins and shapes respectively. Finally, 
Section [6] summarizes the results. 



NclITlG 


Box sizG 


N 


p3,rtr. m&ss 


forcG soft. 


Nhalo 




[Mpc] 






[/i-lkpc] 


> 500 


Wl-20.1 


20 


250^ 


1.36e7 


0.43 


814 


Wl-20.2 


20 


250^ 


1.36e7 


0.43 


849 


Wl-40.1 


40 


250^ 


1.09e8 


0.85 


1132 


Wl-40.2 


40 


2503 


1.09e8 


0.85 


1018 


Wl-90.1 


90 


300^ 


7.19e8 


1.92 


2552 


Wl-90.2 


90 


300^ 


7.19e8 


1.92 


1992 


Wl-90.3 


90 


600^ 


8.98e7 


0.85 


12766 


Wl-180 


180 


300^ 


5.75e9 


3.83 


2406 


Wl-300 


300 


400"^ 


1.12cl0 


4.72 


6415 


W3-20.1 


20 


250^ 


1.32e7 


0.43 


960 


W3-20.2 


20 


250^ 


1.32e7 


0.43 


1011 


W3-40.1 


40 


250^ 


1.05e8 


0.85 


1108 


W3-90.1 


90 


300^ 


6.94e8 


1.92 


1963 


W3-90.2 


90 


300^ 


6.94e8 


1.92 


2043 


W3-90.3 


90 


600^ 


8.67e7 


0.85 


13143 


W3-180 


180 


300^ 


5.55e9 


3.83 


2105 


W3-300 


300 


400^ 


l.OSelO 


4.78 


5255 


W3-360 


360 


400^ 


1.87el0 


5.74 


4947 


W3-360.2 


360 


600^ 


5.55e9 


3.83 


17582 


W5-20.1 


20 


250^ 


1.37e7 


0.43 


974 


W5-40.1 


40 


250^ 


1.09e8 


0.85 


1119 


W5-90.1 


90 


300^ 


7.21e8 


1.92 


1998 


W5-180 


180 


300^ 


5.77e9 


3.83 


2302 


W5-300 


300 


400^ 


1.13el0 


4.74 


5845 



2 N-BODY SIMULATIONS 

Table [2] lists all the simulations used in this paper. Each 
simulation has a unique name, Wx-LL.n, where a; = 1,3, 5 
refers to the WMAPl, WMAP3 and WMAP5 cosmology, 
respectively, LL is the simulation box size in Mpc, and n is 
a sequential integer. For each cosmology we have run simu- 
lations for five different box sizes, which allows us to probe 
halo masses covering the entire range lO^^/i^^ M© ^ M < 
lO^^h-^Mo. In addition, in some case we have run multiple 
(up to three) simulations for the same cosmology and box 
size, in order to test for the impact of cosmic variance (and 
to increase the final number of dark matter haloes). 

All simulations have been performed with PKDGRAV, 
a tree code written by Joachim Stadel and Thomas Quinn 
(Stadel 2001). The code uses spline kernel softening, for 
which the forces become completely Newtonian at 2 soften- 
ing lengths. Individual time steps for each particle are cho- 
sen proportional to the square root of the softening length, e, 
over the acceleration, a: /Sti — r\^JTfai. Throughout, we set 
77 = 0.2, and we keep the value of the softening length con- 
stant in co-moving coordinates during each run. The phys- 
ical values of e at 2 = are listed in Table [2] Forces are 
computed using terms up to hexadecapole order and a node- 
opening angle d which we change from 0.55 initially to 0.7 
sX z = 2. This allows a higher force accuracy when the mass 
distribution is nearly smooth and the relative force errors 
can be large. The initial conditions are generated with the 
GRAFIC2 package (Bertschinger 2001). The starting red- 
shifts Zi are set to the time when the standard deviation 



Cosmological model and Dark Matter halos 3 




10» 10'» 10" 10'2 10'= lO'" 10>s 10" 



M [h-1 Mq ] 

Figure 1. Mass functions for the WMAP3 simulations {only one sim- 
ulation per box size is shown). The doted and solid lines represents the 
Shcth fc Tormen (2001) prediction for the WMAPl and WMAP3 models 
respectively. Error bars indicate the Poisson noise in each mass bin. 



of the smallest density fluctuations resolved within the sim- 
ulation box reaches 0.2 (the smallest scale resolved within 
the initial conditions is defined as twice the intra-particle 
distance). 

In all simulations, dark matter haloes are identified 
using a spherical overdensity (SO) algorithm. Candidate 
groups with a minimum of Nf = 250 particles are selected 
using a FoF algorithm with linking length (p = 0.2 x d (the 
average particle separation). We then: (i) find the point C 
where the gravitational potential is minimum; (ii) determine 
the radius f of a sphere centered on C, where the density 
contrast is A, with respect to the critical density of the Uni- 
verse, pcrit = SH^ /SttG. Using all particles in the corre- 
sponding sphere we iterate the above procedure until we 
converge onto a stable particle set. For each stable particle 
set we obtain the virial radius, rvir, the number of parti- 
cles within the virial radius, A'vir, and the virial mass, Mvir. 
For our adopted cosmologies A ~ 96.7 (WMAPI), A ~ 93.5 
(WMAP3) and A ~ 95.1 (WMAP5). These values are based 
on the fitting function of Mainini et al. (2003). 

Throughout this paper we only use haloes with A^vir > 
500, of which we there are 29944, 50177, and 12238 in the 
combined WMAPl, WMAP3 and WMAP5 simulations, re- 
spectively. For completeness. Fig. [T] shows the halo mass 
function obtained from the WMAP3 simulations, together 
with the predictions from Sheth & Tormen (2001) for the 
WMAPl and WMAP3 cosmologies (error bars show the 
Poisson noise in each mass bin). Note that the mass func- 
tion obtained from our simulation is in excellent agreement 
with this prediction and that the WMAP3 cosmology pre- 
dicts fewer haloes per co- moving volume than the WMAPl 
cosmology, especially at the massive end. 



2.1 Halo parameters 

For each SO halo in our sample we determine a set of param- 
eters, including the virial mass and radius, the concentration 
parameter, the angular momentum, the spin parameter and 
various axis ratios (shape). Below we briefly describe how 
these parameters are defined and determined. A more de- 
tailed discussion can be found in M07. 

2.1.1 Concentration parameter 

To compute the concentration of a halo we first determine 
its density profile. The halo center is defined as the location 
of the most bound halo particle, and we compute the den- 
sity (pi) in 50 spherical shells, spaced equally in logarithmic 
radius. Errors on the density are computed from the Pois- 
son noise due to the finite number of particles in each mass 
shell. The resulting density profile is fit with a NFW profile 
(Navarro et al. 1997): 

pjr) ^ 5c 

Pcrit (r/r,){l + r/r,r' ^' 

During the fitting procedure we treat both rs and 5c as free 
parameters. Their values, and associated uncertainties, are 
obtained via a x minimization procedure using the Lev- 
enberg & Marquart method. We define the r.m.s. of the fit 
as: 

1 ^ 

Prms = (Inpi - Inpm)^ (2) 

i 

where pm is the fitted NFW density distribution. Finally, we 
define the concentration of the halo, Cvir = r^ir/ra, using the 
virial radius obtained from the SO algorithm, and we define 
the error on logc as (ar^/ra) / \n{10), where is the fit- 
ting uncertainty on rg. In order to compare our results with 
previous studies it is also useful to define, C200 = r2oo/rs, 
where r20D is the radius inside of which the density of the 
dark matter halo is 200 times pcrit; accordingly we also de- 
fine M200 and N200 as the mass and the number of particles 
within r2oo. 

2.1.2 Spin parameter 

The spin parameter is a dimensionless measure of the 
amount of rotation of a dark matter halo. We use the defi- 
nition introduced by Bullock et al. (2001b): 



^MvirKir ' vir 

where Mvir, Jvir and Vvir are the mass, total angular momen- 
tum and circular velocity at the virial radius, respectively. 
See M07 for a detailed discussion and for a comparison of 
the different definitions for the spin parameter. 

2.1.3 Shape parameter 

Determining the shape of a three-dimensional distribution 
of particles is a non-trivial task (e.g., Jing & Suto 2002). 
Following AUgood et al. (2006) we determine the shapes of 
our haloes starting from the inertia tensor. As a first step, 
we compute the halo's 3x3 inertia tensor using all the parti- 
cles within the virial radius. Next we diagonalize the inertia 
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tensor and rotate the particle distribution according to tiie 
eigen vectors. In this new frame (in which the moment of in- 
ertia tensor is diagonal) the ratios 3 = 03/01 and p = 03/02 
(where oi 02 ^ 03) are given by: 



03 

Oi 




(4) 



Next we again compute the inertia tensor, but this time 
only using the particles inside the ellipsoid defined by oi, 02 
and 03. When deforming the ellipsoidal volume of the halo, 
we keep the longest axis (oi) equal to the original radius of 
the spherical volume (rvir)- We iterate this procedure until 
we converge to a stable set of axis ratios. 



2.1.^ Offset parameter 

As in M07, for each halo we compute the ofi'set parameter, 
XoS, defined as the distance between the most bound particle 
(used as the center for the density profile) and the center of 
mass of the halo, in units of the virial radius. This ofi'set is a 
measure of the dynamical state of the halo: relaxed haloes in 
equilibrium will have a smooth, radially symmetric density 
distribution, and thus an offset parameter that is virtually 
equal to zero. Unrelaxed haloes, such as those that have only 
recently experienced a major merger, are likely to have a 
strongly asymmetric mass distribution, and thus a relatively 
large XoS ■ Although some unrelaxed haloes may have a small 
XoS, the advantage of this parameter over, for example, the 
actual virial ratio, 2T/V , as a function of radius (Maccio, 
Murante & Bonometto 2003; Shaw et al. 2006), is that the 
former is trivial to evaluate. 
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Figure 4. Histog rams of residuals from the mean C200 rnass relations 
in Fig. [3] Only results for WMAPl and WMAP3 are shown; results for 
WMAP5 are listed in Table [X2I The vertical lines show the 2.3th, 13.9th, 
50th, 84.1th, and 97.7th percentiles of the concentration residuals. The solid 
red line shows a Gauss-Hermite polynomial expansion up to fourth order 
(whose parameters are given in Table rA2t , the dashed red line shows the 
Gaussian corresponding to the zeroth order of this expansion. The logarith- 
mic variance {ctjj-, c) of the Gaussian fit is reported in the top left corner of 
each panel. 



2.2 Relaxed Haloes 

Our halo finder (and halo finders in general) does not dis- 
tinguish between relaxed and unrelaxed haloes. There are 
many reasons why we might want to remove unrelaxed ha- 
los. First and foremost, unrelaxed haloes often have poorly 
defined centers, which makes the determination of a radial 
density profile, and hence of the concentration parameter, an 
ill-defined problem. Moreover unrelaxed haloes often have 
shapes that are not adequately described by an ellipsoid, 
making our shape parameters ill-defined as well. 

One could imagine using prms (the r.m.s. of the NFW 
fit to the density profile) to decide whether a halo is relaxed 
or not. However, while it is true that prms is typically high 
for unrelaxed haloes, haloes with relatively few particles also 
have a high prms (due to Poisson noise) even when they are 
relaxed (cf. Fig. 2 of M07 for the correlation between prms 
and A^vir). Furthermore, since the spherical averaging used 
to compute the density profiles has a smoothing effect, not 
all unrelaxed haloes have a high prms. However, these haloes 
are often characterized by a large offset parameter, x^s- We 
therefore use both prms and XoS to judge whether a halo 
is relaxed or not. Following M07 we split our halo sample 
in unrelaxed and relaxed haloes. The latter are defined as 
the haloes with prms < 0.5 and XoS < 0.07. About 70% of 
the haloes in our sample qualify as relaxed haloes. In what 
follow, we will present results for two different samples of 
haloes: ALL, which includes all haloes with A'vir > 500, 



and RELAXED, which is the corresponding subsample of 
relaxed haloes. 

3 CONCENTRATION-MASS RELATION 

3.1 Differences between WMAPl, WMAP3 and 
WMAP5 

We first discuss the concentration-mass (hereafter c-M) re- 
lation. Figs. [2]and[3]show the Cvir — Afvir and C200 — Af200 rela- 
tions for the WMAPl, WMAP3 and WMAP5 cosmologies, 
and for both the ALL and RELAXED samples, as indicated. 
The symbols show the mean concentrations (in logarithmic 
space) in mass bins of 0.4 dex width. For all three cosmolo- 
gies the c-M relation is well fit by a single power-law (red 
solid line). For the relaxed haloes, these best-fit power-laws 
relations are given by 

logcvir = 1.051 - 0.099 log(Afvir/[lO^^/i"^M0]) (5) 
logc2oo = 0.917- 0.104 log(Af2oo/[10^^/i~'MQ]) (6) 
for the WMAPl cosmology, and 

logcvir = 0.915 - 0.080 log(Afvir/[lO^^/i~^Af0]) (7) 
logc2oo = 0.769 - 0.083 log(Af2oo/[10^^/i"'AfQ]) (8) 
for the WMAP3 cosmology, and 

logCvir = 0.971 - 0.094 log(Afvir/[lO^^/i"^Af0]) (9) 
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Figure 2. c^^^ vs Ai"vir for WMAPI (left), WMAP3 {center) and WMAP5 (right). The upper panels show all haloes with more than 500 particles within 
'"vie while the lower panels show the "relaxed" haloes. The points show the mean concentration {in log c^ir) bins of width 0.4 dex in mass, the error bar 
shows the Poisson error on the mean. The solid lines represent the median concentration in each mass bin, the dashed and dotted lines show the 15.9, 84.1, 2.3 
and 97.7th percentiles of the distribution. The solid (red) line shows a power-law fit to the c-^ij. — A'/^-ir relation: log Cvi^ — zero -|- slope(log M-i^^j- /12h~^ Mq) 
whose parameters are given in the lower left corner of each panel, and in Table FAII 



logcaoo = 0.830 - 0.098 log(A/2oo/[10'^/i-'MQ]) (10) 



for the WMAP5 cosmology. The errors ol these fitting pa- 
rameters are listed in Table lAll Note that these relations 
differ in both the slope and the zero-point. In particular, 
going from WMAPI to WMAP5 to WMAP3, the slopes be- 
comes shallower and the zero-points become smaller. Com- 
paring the WMAPI and WMAP3 cosmologies, which are 
the extremes in terms of the cosmological parameter val- 
ues (see Table [l}, the difference in the mean logCvir is 0.19, 
0.15, 0.11 dex (i.e., a factor of 1.55, 1.41, 1.29 in Cvir) at 
a halo mass of 10^", 10^^ lO"/i"^M0. A similar trend of 
lower normalization and shallower slopes is also seen go- 
ing towards higher redshift for a given cosmology (e.g. Zhao 
et al. 2003a). This supports the notion that the c-M rela- 
tion reflects the assembly histories of dark matter haloes: the 
fact that WMAP3 haloes are less concentrated than their 
counterparts in a WMAPI (or WMAP5) cosmology, simply 
reflects that haloes assemble later in a universe with lower 
Qm and/or lower as. Note that the three cosmologies also 
differ in the spectral index of the matter power spectrum, n, 
which is also responsible for some of the differences in the 
c-M relations. 



Fig. m shows the distributions of the concentration 
residuals, AlogC2oo, relative to the best-fit power-laws 
C20o(M2oo), for the WMAPI and WMAP3 cosmologies. The 
full set of haloes (upper panels) shows a clear skewness to- 
wards low concentrations. This was already pointed out by 
M07, and is even apparent in the simulations of BOl. As 
discussed by M07, this tail of low concentration parameters 
is due to unrelaxed haloes. Removing these haloes results in 
almost Gaussian distributions in A log C200 , as is apparent 
from the lower panels in Fig. U) Table IA2I lists a number of 
parameters for these distributions, including those for the 
WMAP5 cosmology (not shown in Fig. |4|. Note that the 
scatter these distributions is virtually the same for all three 
cosmologies; hence, although the slope and zero-point of the 
c-M relation is clearly cosmology dependent, the scatter is 
not. 

3.2 Comparison of WMAPI results with the 
Millennium Simulation 

Recently, Neto et al. (2007) analysed the c-M relation of 
dark matter haloes in the Millennium simulation (Springel 
et al. 2005). They confirm many of the results published 
previously in M07, namely: 
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Figure 3. Same as Figure^ but for C200 function of -M200 ■ 



(i) The concentration mass relation (at redshift zero) is 
well described by a single power-law with slope ~ —0.10. 

(ii) The inclusion of un-relaxed haloes biases the zero- 
point low, and increases the scatter. 

(iii) The inclusion of un-relaxed haloes is largely respon- 
sible for creating a false correlation between concentration 
parameter and spin parameter. 

The symbols in Fig. [S] show the same mean log C200 as 
function of M200 as shown in Fig. [3] The solid red lines 
in the left-hand panels indicate the best-fit relation ob- 
tained by Neto et al. (2007) from the Millennium Simula- 
tioiS. For both the full sample of haloes and for the sub- 
samples of relaxed haloes, the Neto et al. (2007) results are 
in excellent agreement with our simulation data over the 
range of haloes probed our simulations. Note that, for com- 
parison, the Millennium simulation only covered the range 
10^^/i"^Mq < M200 < 10^^ fe"^M0. This agreement is ex- 
tremely encouraging because the Millennium simulation was 
run with a different A'^-body code (GADGET, rather then 



^ Note that the cosmology adopted for the Millennium simulation 
is very close, but not exactly the same as for our WMAPl cosmol- 
ogy. For the purpose of our discussion, though, these differences 
are completely negligible. 



PKDGRAV), and Neto et al. (2007) used a different halo 
finder. 



Yet, Neto et al. (2007) questioned the resuhs of M07 at 
small masses. They argued that at least 1000 particles are 
needed to reliably determine halo concentrations, and they 
conjecture that the small box size (20 Mpc) of the high- 
est resolution simulation in M07 will be biased by missing 
large scale power. They conclude that some of the differences 
found in M07 between the Eke et al. (2001) model and their 
simulation results at small scales (at a few 10^'^h~^ M©) are 
unreliable. Note, however, that the differences between the 
Eke et al. (2001) model and the M07 simulations are statis- 
tically significant, even at 10^^h~^ Mq (i.e on scales resolved 
with more than 1000 particles). Furthermore, the specu- 
lation by Neto et al. (2007) that the 20 Mpc boxes give 
biased halo concentrations is not justified: as was clearly 
shown in M07, the scatter in the c-M relation is indepen- 
dent of the large scale environment. To further bolster the 
results of M07 at masses below 1 x 10^^ M©, in this pa- 
per we include new simulations of 90 Mpc boxes with 216 
million particles. These simulations have a particle mass of 
9 X lO^/i"^ Mq, which is an order of magnitude higher than 
that of the Millennium simulation. 
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Figure 5. Comparison between the C200 — iVf200 relations in our simulations and those predicted by the toy models of Bullock ct al. (2001a, dotted and 
dashed lines) and NFW (1997) (long-dashed lines). For the Bullock ct al. models we show the value of K that results in the best fit for galaxy mass haloes, 
for the NFW model we use the same value of C for all models. The solid (red) line shows the concentration-mass relations from the Millennium simulations 
as measured by Neto et al. (2007). See text for further details. 



3.3 Comparison with toy models 

A number of studies have presented analytical models for 
the c-Af relation, that have been calibrated against numer- 
ical simulations, but only over a relatively narrow range in 
halo mass. Here we compare our simulation results, which 
cover 5 orders of magnitude in halo mass, to two of the 
most commonly used models: that by Bullock et al. (2001a; 
hereafter the BOl model), and that by Navarro et al. (1997; 
hereafter the NFW model). Both these models have 2 free 
parameters {F and K for BOl, / and C for NFW) that 
have been tuned to reproduce some simulation results. Fig. [5] 
shows the predictions for the c — M relation for these mod- 
els. BOl advocated a model with F = 0.01 and K = 4.0 
which fit the c — M relation of galaxy sized DM haloes 
(10" - lO"/i"^M0). These parameters were obtained for 
a flat ACDM cosmology with Qm = 0.3 and erg — 1.0. A 
slightly lower normalization, K — 3.4, was advocated by 
M07 for a WMAPl cosmology (see Zhao et al. 2003a and 
Kuhlen et al. 2005 for more details). 

This K — 3.4, F = 0.01 model is shown as the long 
dashed line in the upper left-hand panel (ALL sample) of 
Figure [5] It fits the data from M = 10^°/i"^ Mq to M = 
10^'^^"^ M0, but under-predicts the c ~ M relation at high 



masses. Neto et al. (2007) argue that this is a fundamental 
failure of the BOl model. However, as stressed by BOl (in 
both their paper and in their publicly available source code), 
their model is not expected to work on scales where F x 
^^vir ft A/,. A partial solution is to adopt a lower value for 
F; for instance using F — 0.001 and K — 2.7 yields a good 
fit to our simulation data over the range 10^^ Mq M < 
10"h-^ Mq (dotted line). 

By contrast, the NFW model is in reasonable agreement 
with the c-M relation of relaxed haloes at the high mass end. 
This model, with parameters advocated in the original pa- 
per by Navarro et al. (1997) (/ = 0.01, C = 3000), is given 
by the long-dashed lines. This model matches the slope of 
the c — M relation for the highest masses 
However, it dramatically underpredicts the concentrations 
of low mass haloes; in fact, the NFW model predict a slope 
for the c-M relation which is much shallower than what 
we infer from our simulations, causing the discrepancy to 
increase with decreasing halo mass. Neto et al. (2007) gloss 
over these differences by saying that the differences are small 
compared to the scatter in c, i.e. since the scatter in c is a 
factor of 1.4, the discrepancy between the NFW model and 
the simulation results is less than Icr. However, the relevant 
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Figure 6. Comparison between the C200 — A^200 relations in our simulations and those predicted by our revised version of the BOl model. Our ircw model 
is able to simultaneously fit both the low mass and high mass end of the c — M relation. It provides an excellent fit to the simulations results over more than 
4 orders of magnitude in halo mass. 



parameter to compare with is the error on the mean (or me- 
dian) c. Given the numbers in Neto et al. (2007) we obtain 
a Poisson error of 0.004 dex (i.e. less than 1%) for their low- 
est mass bin around 10^^ h^^ M0. We obtain similarly small 
errors with our data for haloes of the same mass. Thus, we 
argue that despite the large amount of scatter in halo con- 
centrations at fixed halo mass, the discrepancy between the 
NFW model and the simulation results arc very significant: 
at M = lO^^h'^ Mq we find the the NFW model underpre- 
dicts the mean concentration by a factor 1.12 (a discrepancy 
of 0.05 dex in logC2oo), which increases to a factor of 1.38 
(0.14 dex in log C200) at M = 10^°/i~^ M©. Of course system- 
atic errors in fitting halo concentrations may well be larger 
than 1%, which would cause them to dominate the error 
budget. However, the fact that Neto et al. (2007) obtain 
results that are virtually identical to us using a different N- 
body simulation code, a different halo finder, and a difi^erent 
fitting procedure, suggests that these effects arc not going 
to mitigate the problem of the NFW model at low masses. 

Both the BOl and NFW models predict that the c-M 
relation is shallower for the WMAP3 cosmology than for 
the WMAPl cosmology, in quantitative agreement with the 
simulation results. However, the problems at high masses for 
the BOl model and at low masses for the NFW models re- 
main. Another problem for these models is that the normal- 
ization of the concentration mass relation in the WMAP3 
cosmology is lower than expected based on the parameters 
K for BOl and C for NFW that are required to match the 
c — M relation in the WMAPl cosmology. Since the BOl 
model does not fit the c-M relation well at high masses, a 



straight forward minimization gives values of K biased 
high. Thus, to determine the best-fit values of K, we con- 
struct a grid of models with K varying at intervals of 0.1, 
and chose the value of K which provides the best fit (by eye) 
over the widest range of halo masses. For F — 0.01 the range 
in mass that is well fitted is M 10^°"^^ft"^ M©, while for 
F = 0.001 the range in mass is M ~ 10"""/i"^ Mq. 

3.4 A revision of the Bullock model 

We now discuss a modification of the BOl model that re- 
tains the agreement for galaxy mass haloes, but improves 
the agreement for group and cluster sized haloes. 

Both the BOl and NFW models are based on the idea 
that the central densities of dark matter haloes reflect the 
mean density of the universe at a time when the central re- 
gion of the halo was accreting matter at a high rate. There- 
fore haloes with central regions that collapse earlier are ex- 
pected to be denser than haloes that collapse later. Thus 
the first step in the BOl (and NFW) model is to assign a 
redshift of collapse, given a halo of mass Myir at a redshift 
of observation, z. 

BOl define the collapse redshift, Zc, as the redshift at 
which the characteristic mass is equal to a fraction F of the 
halo mass at the observation redshift, z, 

M^Zc) = F M^ir{z) , (11) 

Using the spherical collapse formalism, this characteristic 

mass is defined via 

a{M^, Zc) = a{M^, 0)D{zc) = 1.686 (12) 
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where cr(M, z) is the rms overdensity on mass scale M at 
redshift z, and D{z) is the hnear growth rate. 

For comparison, NFW define the collapse redshift as 
the redshift at which, according to the Press & Schechter 
formalism (Press & Schechter 1974; Lacey & Cole 1993), 
half of the virial mass of the halo was first contained in 
progenitors more massive than a fraction / of the final mass. 

The next step in the BOl (and NFW) model is to fink 
the density of the halo at z to that of the mean density of the 
universe at the collapse redshift. First we define the mass of 
the halo via 



M^ir{z) 



,{z)Avir{z)pu{z), 



(13) 



where Avir is the overdensity of the halo relative to the 
mean density of the universe, pu. Alternatively one can set 
AviT{z)pu{z) = A{z)pcTitiz), where A is the overdensity of 
the halo with respect to the critical density of the universe, 
Pcrit. A common choice is to set A = 200, which results in 
the definition of virial mass and radius being independent 
of cosmology (using units of ft = 1). Then the characteristic 
density of the halo at any epoch (as defined by BOl) can be 
written as: 

3Mvir(2) 3 



Ps{z) = 



Avir(2:)pu(2:), 



(14) 



47rrf(2:) 

where c^ir{z) = rvir{z)/rs{z). 

BOl identified the characteristic density of the halo at 
the observation redshift, z, with the mean density of the 
universe at the collapse redshift via 



Pb{z) = K Avir{z)pu{Zc) 



(15) 



where K is a constant to be determined by calibration 
against numerical simulations. Since Pu(zc) = Puiz){l + 
Zc)^/{1 + z)^, the concentration at the observation redshift 
is given by 



■i,{M^i„z) = K{l + z,)/{l + z). 



(16) 



Our modification to the BOl model is to simply assume 
that the characteristic density of the halo ps is independent 
of redshift (i.e. ps(z) = Ps{zc))- Thus 

1/3 



r(Mvir 



K 



Avir(2c) pu{Zc) 



= Kg{zc,z). (17) 



Avir(z) Pu{z) 

As for the original model, the (free) parameter K, which rep- 
resents the concentration of the halo at the collapse redshift 
Zc, has to be calibrated against numerical simulations. The 
function g{zc,z) specifies the growth of the halo concentra- 
tion between the redshifts of collapse and observation. Note 
that this equation is similar to that of the Eke et al. (2001) 
model, except that they implicitly assumed that Cvir{zc) = 1 
(which is wrong), and they define Zc differently. 

We can also express g in terms of A (2:) and pcrit 

1/3 

(18) 



9(zc, 



A( 



Pcrit ( 2c) 



A(2) p, 

In the case of the A 



rit(2) 



200 halo definition, we thus have that 



C200 (2) = K200 



pcrit (2c) 
Pcrit (2) 



1/3 



H{z,) 
H{z) 



2/3 



(19) 



So that the evolution in C200 is given by the evolution of the 
Hubble parameter, which is given by: 
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log M200 [h-'Mg] 

Figure T. Growth factor of concentration parameter as a function of lialo 
mass (computed at 2 — 0) for WMAPl cosmology. Tire dashed line shows 
the growth factor for the BOl model, while the solid line shows the growth 
factor for our revised BOl model. Our model has a shallower mass depen- 
dence at high masses 10^^ h "'"Mq), resulting in a better agreement 
with the high mass end of the concentration mass relation (cf. Fig, fsl . 



H\z 



"0 



^ = [S^A + n„(l + zf + (1 - Ha - !^m)(l + zf] , .(20) 



Our modification of the BOl model only consists of tak- 
ing account of the redshift dependence of the halo density 
contrast. In the case of the A — 200 halo definition, the red- 
shift dependence of Avir ~ A/Slm is completely governed 
by that of Qm- For an EdS universe Qrn{z) = 1, so that our 
revised model yields exactly the same Cvir(2) as the orig- 
inal BOl model. However, for a ACDM universe Qrn is a 
function of redshift, so that the difference between i}rn{zc) 
and Q,,n{z) can be significant, and thus also the difference 
in Cvir(2) between the BOl model and our revised version. 

FiglZ] compares the concentration growth factors, g{zc) 
for the BOl model {g{zc) = (1 -I- 2c)) and for our new model 
{g{zc) = [H{zc)/Hof''^) in a WMAPl cosmology. For haloes 
below M = lO^^/i"^ M0 (and collapse redshifts > 2) the 
slopes of g{M2oo) are the same for both models. This owes 
to the fact that Q,rn{z) — * 1 at high redshifts. However for 
massive haloes (and low collapse redshifts), the concentra- 
tion growth factor g has a shallower slope in our model than 
in the original BOl formulation. Because of this difference 
our new model predicts concentrations at the high mass end 
that are relatively higher than for the original BOl model. 
Note also that at the low mass end, the new concentration 
growth factor is lower than for the BOl model; consequently, 
in order to match the same c-M relation at the low mass end, 
the K parameter in our new formulation has to be somewhat 
higher than in the BOl model. 

The symbols in Fig. [6] show the c-M relation obtained 
from our simulations for all three cosmologies, separately for 
all haloes (left-hand panel) and for the subsample of relaxed 
haloes (right-hand panel) . The lines correspond to the best- 
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Figure 8. Formation rcdshifts as a function of lialo mass for tire BOl and 
NFW models. The NFW model has a much weaker mass dcpeirdenee to the 
formation redshift than the BOl model, which results in a c — A/ relation 
much shallower than in the simulations as shown in Figure fsl 



fit c-M relations obtained using our new toy model, in which 
we let K hea, free parameter. The best-fit values of K are in- 
dicated. Note that our new model can fit the simulation data 
remarkably well over the entire range of halo masses probed, 
and for all three cosmologies. Clearly, this new model is a 
significant improvement over the BOl and NFW models. As 
with the NFW model and the original BOl model, however, 
different cosmologies require a different normalization pa- 
rameter K. In particular, the best-fit value of K decreases 
going from WMAPl to WMAP5 to WMAP3. This is un- 
fortunate, as there is currently no existing model that can 
a priori predict the value of K given the cosmology. Conse- 
quently, for each cosmology, the best-fit value of K has to 
be determined empirically using high-resolution numerical 
simulations. For the three cosmologies considered here, the 
differences in the best-fit values of K are small (~ 5%), but 
they nevertheless indicate that our model does not fully cap- 
ture the power spectrum dependence of the c — M relation 
(nor does the original BOl model. We plan to address this 
issue in more detail in a future paper. 



3.5 Why does the NFW model fail? 

There are two differences between the NFW model on the 
one hand and the BOl model and it modification suggested 
here on the other: the definition of halo formation redshift, 
and the method used to link the mean density of the universe 
at the formation redshift to the concentration of the halo at 
the observation redshift. Fig,|S]shows the formation redshifts 
as a function of halo mass for the NFW and BOl models in 
a WMAPl cosmology. The NFW model predicts formation 
redshifts with a significantly shallower mass dependence that 
the BOl model. We have experimented with a model using 
the NFW formation time, but the BOl method of linking 



formation redshift to halo concentration. For none of these 
models, however, were we able to obtain a concentration- 
mass relation with a slope as steep as found in the simu- 
lations. Consequently, we conclude that the main problem 
with the NFW model is its definition of the halo collapse 
redshift. 



3.6 Comparison with observations 

Several studies have addressed the issue of the apparent in- 
consistency between the inner density slopes of dwarf and 
LSB galaxies and those predicted by CDM (e.g., Moore 
1994; Flores & Primack 1994; van den Bosch & Swaters 
2001; de Blok et al. 2001, 2002; Swaters et al. 2003, Button 
et al. 2005; Gentile et al. 2005; Simon et al. 2005; Kuzio de 
Naray et al. 2008). For more massive spiral galaxies the sit- 
uation is more complicated because the baryons contribute 
significantly (and may even dominate) the total mass pro- 
file in the inner regions, and thus pure CDM predictions are 
difficult to be tested. For these reasons we limit our com- 
parison to data of dwarf and LSB galaxies. Measuring inner 
density slopes of dark matter haloes from rotation curves is 
difficult. A common practice is to fit a power-law to the rota- 
tion curve, either to the inner few data points (e.g., de Blok 
2001); or to the full rotation curve (e.g. Simon et al. 2005). 
Both methods have some drawbacks. First of all, the in- 
ner few points of the rotation curve carry with them the 
largest systematic errors. In addition, realistic dark matter 
halo density profiles (whether they be pseudo-isothermal, 
NFW, or generalized NFW profiles) all have continuously 
varying density slopes with radius, so that the results of 
fitting a single power-law are extremely sensitive to the ex- 
act radial range covered. A more robust approach is to fit a 
double power-law density profile to the entire rotation curve 
(e.g., van den Bosch & Swaters 2001, Button et al. 2005). 
However, as emphasized by these studies, even in this case 
there are certain degeneracies inherent to the mass mod- 
eling that make it difficult to obtain stringent constraints 
on the inner density slopes, even with data of high spatial 
resolution. 

In addition to predictions for the slope of the inner 
density profile, the CBM model also makes predictions re- 
garding the concentration of dark matter haloes. Unfortu- 
nately, measuring halo concentrations requires knowledge of 
the virial radius of the halo, a parameter which is poorly 
constrained from the observations. This is due to the fact 
that a rotation curve only probes the inner ~ 10% of the 
dark matter halo. 

An observationally more robust measurement of the 
central density is given by the dimensionless parameter 
Ay/2, defined as the average density of the halo, with re- 
spect to the critical density, inside the radius Rv/2, where 
the halo circular velocity drops to half its maximum value 
(Alam, Bullock & Weinberg 2002). 



Av/2 = 



V/2) 



Pcrit 



= 50 



Vmax 


2 


h ^kpc 


. kms~^ . 




Rv/2 



(21) 



For an NFW halo, the circular velocity reaches a maxi- 
mum, Vmax, at a radius r ~ 2.163 r^, so that Rv/2 ~ 0.126rs. 
The maximum circular velocity is given by 

VL^ = 0.2162 Kl Cvir//(Cvir), (22) 
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Figure 9. Comparison between the A^^2 ~ ^max relation from obser- 
vations of dwarf and LSB galaxy rotation curves and the predictions of 
N-body simulations in a variety of ACDM cosmologies. The observations 
are: de Blok, McGaugh & Rubin (2001, dBOl); de Blok & Bosma (2002, 
dB02); and Swaters et al. (2003, S03). The shaded region shows the 68.3% 
and 95.4% range of A^y2 from our WMAP5 simulations. See text for fur- 
ther details. 



where f{x) — ln(l -|- a;) — a::/(l -I- x). Note that the factor 
0.2162 is equal to f{x)/x, with x = 2.163. The virial velocity 
scales with the virial mass and virial radius according to 



Kir 



kms" 



Kir 

kms" 



G 



1/3 



Mvir 


1/3 


r A 1 


1/6 






[200 J 





h ^kpc 



200 



1/2 



(23) 



(24) 



with G the gravitational constant. Using these relations 
and the definition of the concentration parameter (cvir = 
'"vir/fs), can be expressed purely in terms of the con- 

centration parameter: 



^V/2 



680.9 



A 



200 /(Cvir 



(25) 



Thus given a relation between Cvir and Mvir for NFW haloes, 
we can convert this into a relation between Ay/2 and Vmax- 



Fig. 1^ shows the relation between Ay/2 Knax for 
both observations (symbols) and theory (lines). The obser- 
vations are for dwarf and LSB galaxies, color coded accord- 
ing to the reference. These values are calculated based on the 
pseudo-isothermal halo fits to the observed rotation curves. 
We have removed a few galaxies with obviously very bad fits, 
or for which Kiax was poorly constrained by the data. We 
have verified that including these galaxies makes no signifi- 
cant difference to the median observed Ay/2. The solid black 
circle shows the median of the data points, and the error bars 
reflect the median errors on Ay/2 and Kiax. Although these 
galaxies are dark matter dominated, the baryons make a 



circle shows the median A 



V/2 



and Knax for the dark halo 



when the contribution of the baryons (stars and gas) to the 
rotation curve has been subtracted. These baryonic contri- 
butions are available for about half of the galaxies. This re- 
sults in a Ay/2 that is approximately 0.16 dex (i.e., a factor 
of 1.45) lower, and a Kiax that is ~ 0.04 dex (i.e., factor of 
1.10) lower. We note that the median Ay/2 and Knax with- 
out subtracting the baryons is the same for this subset, as 
the full sample, so that the differences in Ay/2 and Knax are 
not simply caused by a selection effect. 



The lines show the predictions for A 



V/2 



Kiax for 4 



cosmologies. The dotted line shows the ACDM cosmology 
used in Alam et al. (2002), which has fim = 0.3, Q.k ~ 0.7, 
0-8 = 1.0, h = 0.7 and n = 1. As noted by Alam et al. (2002), 
this cosmology results in central densities that are a factor of 
4 higher than observed. However, this cosmology has values 
of (jg, fim and n that are significantly higher than those for 
the WMAP cosmologies considered here. The Ay/2 — Knax 
relations for WMAPl, WMAP3 and WMAP5 cosmologies 
are shown as the short dashed, long-dashed and solid lines, 
respectively. Note that the data are broadly consistent with 
all these cosmologies; within the errors, the data seem to 
prefer a cosmology with erg ~ 0.8 and n ~ 0.96. 

Note, though, that there are several systematic effects 
that could alter this conclusion. If the effect of halo contrac- 
tion (e.g. Blumenthal et al. 1986) is taken into account, it 
will result in a lower Ay/2 for the initial halo and, and thus 
favor cosmologies with lower ag. On the other hand, there 
are several systematic effects that result in underestimates 
of Kiax (see Swaters et al. 2003; Rhee et al. 2004; Ifayashi 
& Navarro 2006; Valenzuela et al. 2007), which will cause 
the opposite effect and thus would favor cosmologies with 
higher erg. 

For the WMAP5 cosmology we show the 1 and 2 — cr 

- M rela- 
; 0.28 dex 
is ~ 0.43 



dispersion in Ay/2 assuming a scatter in the c 
tion of (Jinc = 0.26 (which results in a scatter of : 

v/2 



in Ay/2). The dispersion of the observed A 



dex. Ifowever, subtracting a measurement uncertainty of 
0.32 dex (the median measurement uncertainty from the ob- 
servations) results in an intrinsic scatter of ~ 0.29 dex, in 
good agreement with our theoretical predictions. 

We conclude therefore, that modulo the caveats of halo 
contraction and systematic effects (which to first order tend 
to cancel each other), the central densities of dwarf and LSB 
galaxies are in good agreement with predictions for a ACDM 
cosmology with parameters favored by the WMAP mission. 

However, on larger scales observations seem to require 
higher concentrations than predicted by ACDM . Using X- 
ray ( Chandra and XMM-Newton) observations of 39 massive 
galaxy clusters, Buote et al. (2007) obtained Cvir = 9.0 ±0.4 
at Mvir = 10^'* M© for their full sample (covering the mass 
range lO" < Mvir < IO^'^Mq), and Cvir = 7.6 ± 0.5 for the 
sample restricted to Mvir 5S lO^^M©. Combining strong lens- 
ing and X-ray observations, Comerford & Natarajan (2007) 
found Cvir = 11.1 at Afvir = lO^'*Af0. In our simulations 
we find the median concentrations of relaxed haloes with 
Afvir = 10"Mq to be Cvir = 6.9 (WMAPl) and Cvir = 5.5 
(WMAP3). 

The median Cvir in the full sample of Buote et al. is 
a factor 1.30 higher than in our WMAPl simulations and 
a factor 1.63 higher than in our WMAP3 simulations. For 
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Figure 10. Spin parameter vs. mass for all halos {upper panels) and 
relaxed halos (lower panels). WMAP5 result arc not shown since they do 
not present any significant difference from WMAPl or WMAP3. The points 
represent the median spin in each mass bin, the error bar shows the Poisson 
error on the mean. The dashed and dotted lines show the 15.9, 84.1, 2.3 and 
97.7th percentiles of the distribution. The solid (red) line shows a power- 
law fit to the A — M2QQ relation; log A — zero + slope(log M2QQ /^'2h~^ Mq) 
whose parameters are given in the lower left corner of each panel. Linear 
fit parameters for the WMAP5 model are reported in Table FAII 



Figure XI. Histograms of the distribution of the halo spin parameters. 
Only results for WMAPl and WMAP3 are shown, results for WMAP5 are 
listed in Table 1X2] The vertical lines show the 2.3th, 13.9th, 50th, 84.1th, 
and 97.7th percentiles of the spin residuals. The solid red line shows a 
Gauss-Hermite polynomial expansion up to fourth order (whose parameters 
are listed in Table lA2t . the dashed red line shows the Gaussian correspond- 
ing to the zeroth order of this expansion. The mean and the logarithmic 
variance (fT^^ a) °f Gaussian fit are reported in the top left corner of 
each panel. 



the high mass sample of Buote et al. the discrepancies are 
smaller: a factor of 1.1 for WMAPl and 1.25 for WMAP3. 
Given that the scatter in halo concentrations is a factor of 
1.3, selection biases in the data (such as selecting the most 
relaxed clusters, which are likely to form earlier and thus to 
have higher concentration) could plausibly reconcile the ob- 
served concentrations from Buote et al. with the WMAP3 
cosmology. The concentrations of Comerford & Natarajan 
(2007) are a factor 1.6 higher than our WMAPl results and 
a factor 2.0 higher than WMAP3 results, and selection ef- 
fects would have to be rather severe to reconcile even our 
WMAPl results with these observations. Another possibility 
is that halo contraction due to the condensation of baryons 
at the center of the halo (Blumenthal et al. 1986) has played 
an important role. However, studies which model the ra- 
dial density profiles of clusters and elliptical galaxies suggest 
that it is difficult to reconcile models with adiabatic contrac- 
tion with observations (Zappacosta et al. 2006; Humphrey 
et al. 2006, Gastaldello et al. 2007). Thus, it seems that there 
is a discrepancy between model and data at the high mass 
end, but not at the low mass end. 

This may signal the need for different normalizations of 
the power spectrum on different scales, which in turn may 
indicate that the power spectrum has a significant tilt, or 
a running spectral index with a shallower slope (lower n) 
at lower masses. However, systematic errors and selection 
effects in the data need to be understood better before such 
a conclusion can be verified. 



4 SPIN PARAMETERS 

Fig. [To] shows the spin parameter versus halo mass for 
the WMAPl (left) and WMAP3 (right) cosmologies for all 
haloes (top) and for the subsample of relaxed haloes (bot- 
tom) . As in M07, we find no mass or cosmology dependence 
of the halo spin parameters. For this reason we do not in- 
clude the results for the WMAP5 cosmology in the plots. 
For completeness, we do list the corresponding parameters 
in Tables lAll and lA2l in the Appendix. 

Fig. [TT] shows the distribution of spin parameters. As 
with the scatter in concentration, the scatter in A is not 
perfectly log-normal. The distribution of spins for all haloes 
shows a small skewness to low spin parameters. However, we 
caution against an over interpretation of this skewness, as 
Bullock et al. (2001b) have demonstrated that haloes with 
lower spins have larger associated uncertainties, with the 
error in A given roughly by 0.2A"^iV"^''^. Thus a halo with 
1000 particles will have an uncertainty of 63%, 18%, and 6% 
for a spin parameter of 0.01, 0.035, and 0.10. Even for haloes 
with 10000 particles the error on a spin parameter of 0.01 is 
20%, which will introduce a non-negligible skewness to the 
distribution of A. 

We find that the mean, dispersion and skewness of the 
distributions (see Table IA2|) are remarkably similar for all 
three cosmologies considered here, both for the full set of 
haloes and for the subsample of relaxed haloes. The relaxed 
haloes have a smaller mean and variance, and a larger skew- 
ness. This is due to the removal of un-relaxed haloes which 
typically have higher spins than relaxed haloes (see M07 for 
details). 
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Figure 12. Short to long axis ratio, 3, vs mass for WMAPl (left) and 
WMAP3 (right), and for all haloes (upper) and relaxed haloes (lower) with 
-^200 > 3000. WMAP5 model results (not shown) lie in between WMAPl 
and WMAP3. The points represent the median value of s spin in each mass 
bin, the error bar shows the Poisson error on the mean. The dashed and 
dotted lines show the 15.9, 84.1, 2.3 and 97.7th percentiles of the distribu- 
tion. The solid (red) line shows a power-law fit to the s — A/200 relation: 
s — zero + slope(log A-/200/ A/q ) whose parameters are given in the 
lower left corner of each panel. Linear fit parameters for the WMAP5 model 
are reported in Tabic IaTI 



5 HALO SHAPES 

Fig [12] shows the relation between s (defined as the ratio 
between the short and long axes) and M2oo- Here we only 
consider haloes with at least 3000 particles, because with 
fewer particles we find evidence for resolution effects. The 
s — M200 relation is well fitted with a power-law of slope 
~ —0.05, in all cosmologies and for all and relaxed haloes. 
Fitting parameters for the shape-mass relation and for its 
scatter, for all three cosmological models, are listed in Tables 
[KT\ and 1X2] in the Appendix. For the WMAP3 cosmology, 
the zero point for relaxed haloes is 0.03 higher than for the 
full sample of haloes. In addition, we find that the zero-point 
if 0.05 higher for the WMAPl cosmology compared to that 
of the WMAPl cosmology. On the other hand, the relation 
between p (defined as the ratio between the short and inter- 
mediate axes) and A/200 shows a much weaker correlation, 
and only a marginal dependence on cosmology (Fig. I13p . 

Figs. [14] & [15] show the distributions of s and p about 
the mean relations shown in Figs. [12] & [13] These are 
roughly Gaussian. The full sample reveals a mild skewness to 



Figure 13. Same of Fig.[T2]but for the middle axis, p, vs mass. WMAP5 
model results {not shown) lie in between WMAPl and WMAP3. Values for 
the slope and the zero of the linear fit for the WMAP5 model are reported 
in Table IaTI 




Figure 14. Histo grams of scatter in the halo shape parameter s (ratio 
between minor and major axis length) for WMAPl (left) and WMAP3 
(right) ; results for WMAP5 (not shown in the plot) are listed in Table 
[A2] The upper panels show all haloes with ^^200 > 3000, while the lower 
panels show the relaxed haloes. The dashed (red) lines show the best fitting 
Gaussian, while the solid (red) lines show the Gauss-Hermite polynomial 
expansion (whose parameter are reported in Table [A2 t . The variance of the 
Gaussian fit is reported in the upper left corner of each panel. 
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Figure 15. Histograms of scatter in the halo shape parameter p (ratio be- 
tween minor and intermediate axis length) for WMAPl (left) and WMAP3 
(right); results for WMAP5 (not shown in the plot) are listed in Table [A2] 
The upper panels show all haloes with greater than 3000 particles within 
^200' while the lower panels show the "relaxed" haloes. The dashed (red) 
lines show the best fitting Gaussian, while the solid (red) lines show the 
Gauss-Hcrmite polynomial expansion (parameters listed in Tablc[A2}. The 
variance of tlie Gaussian fit is reported in tire upper left corner of each 
panel. 



Figure X6. Shape within O-Sr^ir versus mass in different eosmological 
models. A universe with lower erg and Vljn produces haloes that are more 
elongated. In the upper panels the dotted lines give the relationship sug- 
gested by Allgood et al. (2006). The solid lines show a fit of the form shape 
= zero + slope(log M^i,/M, ), where log M» = 12.82, 12.17/i~^ Mq in the 
WMAPl and WMAP3 cosmologies, respectively. The slopes and zero points 
of these fits are given in the lower left corner of each panel. The long-dashed 
lines show the corresponding relations for the shapes measured within the 
virial radius. The short dashed lines enclose 68.3% of the points, and cor- 
respond to a scatter of ^ 0.1. 



low s, which is most likely due to the presence of unrelaxed 
haloes. 

Allgood et al. (2006) presented a detailed analysis of 
the dependence of halo shape on mass and on the underly- 
ing eosmological model. They found that the redshift, mass 
and erg dependence of the mean smallest-to-largest axis ra- 
tio of haloes is well described by a simple power law relation 
(so.a) = a{Mvii/Mt)^ , where in this case so. 3 is measured in- 
side O.Srvir, and the z and as dependencies are governed by 
the characteristic non-linear mass, M* = M,{z,as). Using 
several simulations they found the following values for the 
fitting parameters: a = 0.54 ± 0.02 and b = -0.050 ± 0.003. 



Fig. [TC] shows the mean so.3-mass and po.3-niass rela- 
tions in our simulations. The dotted line shows the relation 
found by Allgood et al. (2006) where we used the following 
values for the characteristic mass: log(M,/h~^M0) =12.82, 
12.17, for the WMAPl and WMAP3 cosmologies, respec- 
tively. Our simulations show a similar slope, but somewhat 
higher normalization compared to Allgood et al. (2006). In 
addition, we find that the differences in the shapes between 
WMAPl and WMAP3 cosmologies are not simply explained 
by a difference in A/* , as proposed by Allgood et al. (2006) . 
This can be seen in the upper panels of Fig.fTHlby comparing 
the difference between the Allgood et al. (2006) prediction 
(dashed line) with the results of our simulations (solid line). 
This difi'erence is not constant and increases from WMAPl 
to WMAP3. 

The dot-dashed lines in Fig [16] show s and p within the 



virial radius. This shows that s decreases towards the center 
of the halo, in good agreement with Allgood et al. (2006). 
The dispersion in the axis ratios measured at O.Srvir is 
0.10, only marginally smaller than the dispersions measured 
within Tvir. Results for the shape-mass relation and for the 
shape distribution for the WMAP5 model (not shown in the 
figures) are summarized in Table \KT\ and in Table IX2I 



6 SUMMARY 

In this paper we have used a large set of eosmological N-body 
simulations in WMAPl, WMAP3 and WMAP5 cosmologies 
to study how changes in the eosmological parameters affect 
the halo mass dependence of the concentration parameter, c, 
the spin parameter. A, and the halo shape parameters, s and 
p. The simulations span 5 orders of magnitude in halo mass 
(10^" — W^^h~^ M0), covering the entire range from haloes 
those that host individual dwarf galaxies to those associated 
with massive clusters. 

At a fixed mass, haloes in a WMAP3 cosmology are 
significantly less concentrated than their counterparts in a 
WMAPl cosmology. As already noted by other authors this 
can be ascribed to the lower value for us in the WMAP3 cos- 
mology compared to WMAPl, that shifts the entire process 
of structure formation towards lower redshifts. In particu- 
lar, for a halo of lO^^ft'^M© the WMAP3 concentrations 
are 1.41 times lower than in a WMAPl cosmology. Due to 
this lower normalization, the central densities of dark mat- 
ter haloes in the WMAP3 cosmology are consistent with the 
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observed central densities of dark matter haloes of dwarf 
and LSB galaxies. However, on the scale of clusters, the 
WMAP3 concentrations may actually be too low compared 
to observational constraints from X-ray measurements and 
from gravitational lensing. The WMAP5 concentrations are 
intermediate between the WMAPl and WMAP3 cosmolo- 
gies: they are still consistent with the data on dwarf and 
LSB galaxies, but somewhat too low on the scale of galaxy 
clusters. Although this may indicate a problem with the ex- 
act shape of the power spectrum of density perturbations, 
systematic errors and selection effects in the cluster data 
need to be better understood, before the data can be used 
to constrain the cosmological spectral index. 

We find that for all three cosmologies, the average halo 
concentration as function of halo mass is well fitted by a 
single power-law over the entire range of halo masses cov- 
ered by our simulations. This is inconsistent with all exist- 
ing models for the mass dependence of halo concentrations, 
which predict that the slope of the c-M relation becomes 
steeper at higher masses. In particular, the model suggested 
by Navarro, Frenk, & White (1997) matches the slopes and 
(to a lesser degree) normalizations of the c-M relation at the 
massive end {M > lQ^'^h~^ M©), but dramatically underpre- 
dicts the concentrations for low mass haloes. The model sug- 
gested by Bullock et al. (2001a), on the other hand, matches 
the slopes of the c-M relations for M < lO^^/i"^ M©, but 
underpredicts the concentrations for more massive haloes. 

We propose a modification to the BOl model, based on 
the idea that the characteristic density of a halo remains con- 
stant after the halo forms. This results in a growth in halo 
concentration proportional to H{zc)'^^^, rather than (1-|-Zc) 
as in the BOl model. Consequently, our new model repro- 
duces the slope of the c — M relation over the full range 
of masses in our simulations, and for each cosmology. How- 
ever, as for the original BOl model, the normalization is 
cosmology dependent (at the level of a few percent), sug- 
gesting that the model still does not completely capture the 
relevant dependencies on the power spectrum. Until such a 
model is available, the normalization of the c-M relations 
will have to be calibrated against numerical simulations for 
each cosmology, if precision concentration parameters are 
required. 

In agreement with previous studies (e.g. Bullock 
et al. 2001b) we find that the distribution of spin param- 
eters is independent of halo mass and cosmology. There is 
a trend, though, that less relaxed haloes have higher spin 
parameters, as previously noticed by Maccio et al. (2007). 
Finally, we find haloes to be more flattened in the WMAP3 
cosmology than in the WMAPl cosmology, consistent with 
the suggestion that haloes that form later are more aspher- 
ical (Allgood et al. 2006). 
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APPENDIX A: PARAMETERS 

In the appendix we summarize all the parameters of the fit- 
ting functions used through this paper. Table lAll contains 
the slope, zero and relative errors for the power-laws fits 
shown in Figs. El 151 [TOI [T^ Table IMl shows the parameters 
describing the distribution around the mean of concentra- 
tions, spin and shapes in the three different cosmological 
models (Figs: H] [TT] [TH [T5)) . While this paper was ready 
for submission a similar study of the concentration mass- 
relation for the WMAP5 cosmology was presented in Duffy 
et al. (2008). They limited their studied only to haloes de- 
fined with more than 10* particles inside the virial radius. 
This left them with only 1269 haloes, nevertheless the slope 
of their power law fit is consistent with our results (within 
the errors) for both the ALL and RELAXED sample. 
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Table Al. Parameters of the power-law fits shown in Figs[2l l3lll0lll2l fc ll3l The data are fitted with the power law y =zero-|-slopc(log M— 
12), where M is in units of h~^M0. The second column reports the mass range over which the fit has been obtained, the third the 
number of haloes within this mass range (using all the simulations listed in Table [2)l- 



sample 




logj^o ^'^ 




zero 


error 


slope 


error 


















log Cvir log A/vir 
















Wl-ALL 


500 


9.83-14.95 


29693 


1.011 


0.001 


-0.114 


0.001 


Wl-RELAXED 


500 


9.83-14.95 


21876 


1.051 


0.001 


-0.099 


0.001 


W3-ALL 


500 


9.82-14.95 


49830 


0.861 


0.001 


-0.086 


0.001 


W3-RELAXED 


500 


9.82-14.95 


33913 


0.915 


0.001 


-0.080 


0.001 


W5-ALL 


500 


9.84-14.86 


12184 


0.925 


0.001 


-0.108 


0.001 


W5-RELAXED 


500 


9.84-14.86 


8282 


0.971 


0.001 


-0.094 


0.001 


log C200 log -A/200 
















Wl-ALL 


500 


9.83-14.92 


25952 


0.879 


0.001 


-0.119 


0.001 


Wl-RELAXED 


500 


9.83-14.85 


19528 


0.917 


0.001 


-0.104 


0.001 


W3-ALL 


500 


9.82-14.82 


40027 


0.719 


0.001 


-0.088 


0.001 


W3-RELAXED 


500 


9.82-14.82 


28344 


0.769 


0.001 


-0.083 


0.001 


W5-ALL 


500 


9.84-14.93 


9988 


0.787 


0.001 


-0.110 


0.001 


W5-RELAXED 


500 


9.84-14.93 


7060 


0.830 


0.001 


-0.098 


0.001 


log A vs log A/200 
















Wl-ALL 


500 


9.83-14.92 


25733 


-1.474 


0.002 


0.005 


0.002 


Wl-RELAXED 


500 


9.83-14.92 


19377 


-1.513 


0.002 


-0.007 


0.002 


W3-ALL 


500 


9.82-14.82 


39662 


-1.474 


0.001 


-0.002 


0.001 


W3-RELAXED 


500 


9.82-14.82 


28114 


-1.526 


0.001 


-0.008 


0.001 


W5-ALL 


500 


9.84-14.93 


9988 


-1.458 


0.004 


0.001 


0.003 


W5-RELAXED 


500 


9.84-14.93 


7060 


-1.505 


0.004 


-0.009 


0.004 


s vs log Af200 
















Wl-ALL 


3000 


10.61-14.92 


4886 


0.654 


0.015 


-0.057 


0.016 


Wl-RELAXED 


3000 


10.61-14.85 


3820 


0.686 


0.017 


-0.056 


0.018 


W3-ALL 


3000 


10.60-14.82 


6022 


0.598 


0.016 


-0.043 


0.013 


W3-RELAXED 


3000 


10.60-14.82 


4273 


0.630 


0.018 


-0.046 


0.016 


W5-ALL 


3000 


10.60-14.82 


1492 


0.623 


0.035 


-0.052 


0.025 


W5-RELAXED 


3000 


10.61-14.93 


1060 


0.657 


0.040 


-0.054 


0.029 


pvs logAf200 
















Wl-ALL 


3000 


10.61-14.92 


4886 


0.820 


0.015 


-0.022 


0.016 


Wl-RELAXED 


3000 


10.61-14.85 


3820 


0.824 


0.017 


-0.018 


0.018 


W3-ALL 


3000 


10.60-14.82 


6022 


0.794 


0.016 


-0.011 


0.013 


W3-RELAXED 


3000 


10.60-14.82 


4273 


0.801 


0.018 


-0.009 


0.016 


W5-ALL 


3000 


10.60-14.82 


1492 


0.807 


0.035 


-0.016 


0.025 


W5-RELAXED 


3000 


10.61-14.93 


1060 


0.812 


0.040 


-0.014 


0.029 
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Table A2. Parameters of the distribution of concentrations (Fig|4j , spin (Fig.[TTJ, and shapes (Figs. [Ml fc llSj l. The distribution have 
been fitted by a Gauss-Hermite polynomial expansion up to fourth order {h\ — hi). The zeroth order of this expansion is a Gaussian fit, 
whose mean < Alogc > and dispersion crAiogc are determined by setting hi = h2 = 0. Columns 6-10 show the 2.3, 15.9, 50.0, 84.1, and 
97.7th percentiles of the distribution. 



sample 


< A log c> 


log c 


h3 


hi 


2.3th 


15.9th 


50th 


84.1th 


97.7th 


A log C200|Af200 




















Wl-ALL 


0.033 


0.129 


-0.670 


0.356 


-0.405 


-0.143 


0.022 


0.146 


0.269 


Wl-RELAXED 


0.011 


0.111 


-0.421 


0.216 


-0.263 


-0.113 


0.007 


0.115 


0.228 


W3-ALL 


0.040 


0.132 


-0.880 


0.252 


-0.417 


-0.147 


0.027 


0.150 


0.252 


W3-RELAXED 


0.015 


0.109 


-0.622 


0.146 


-0.260 


-0.113 


0.010 


0.114 


0.209 


W5-ALL 


0.041 


0.130 


-0.929 


0.351 


-0.417 


-0.148 


0.028 


0.149 


0.252 


W5-RELAXED 


0.015 


0.105 


-0.620 


0.259 


-0.268 


-0.110 


0.011 


0.112 


0.209 


log A 




















Wl-ALL 


-1.470 


0.258 


-0.158 


0.050 


-2.069 


-1.752 


-1.478 


-1.225 


-0.984 


Wl-RELAXED 


-1.499 


0.236 


-0.307 


0.034 


-2.083 


-1.774 


-1.508 


-1.284 


-1.085 


W3-ALL 


-1.465 


0.256 


-0.194 


-0.003 


-2.062 


-1.747 


-1.472 


-1.227 


-1.004 


W3-RELAXED 


-1.514 


0.235 


-0.297 


-0.019 


-2.095 


-1.783 


-1.522 


-1.303 


-1.110 


W5-ALL 


-1.466 


0.253 


-0.162 


0.051 


-2.070 


-1.743 


-1.468 


-1.218 


-0.942 


W5-RELAXED 


-1.508 


0.228 


-0.341 


0.041 


-2.105 


-1.774 


-1.515 


-1.297 


-1.065 


As 1 A/200 




















Wl-ALL 


0.017 


0.121 


-0.343 


0.098 


-0.317 


-0.121 


0.011 


0.130 


0.230 


Wl-RELAXED 


0.004 


0.112 


-0.123 


-0.123 


-0.224 


-0.110 


0.003 


0.112 


0.205 


W3-ALL 


0.007 


0.123 


0.215 


0.169 


-0.269 


-0.121 


-0.001 


0.129 


0.255 


W3-RELAXED 


-0.006 


0.116 


0.441 


0.011 


-0.218 


-0.113 


-0.004 


0.119 


0.232 


W5-ALL 


0.008 


0.128 


-0.117 


0.196 


-0.294 


-0.129 


0.004 


0.133 


0.248 


W5-RELAXED 


-0.002 


0.117 


0.321 


-0.031 


-0.233 


-0.114 


-0.003 


0.123 


0.221 


Ap| A/200 




















Wl-ALL 


0.011 


0.103 


-0.779 


-0.514 


-0.209 


-0.102 


0.006 


0.103 


0.162 


Wl-RELAXED 


0.011 


0.100 


-0.867 


-0.504 


-0.203 


-0.099 


0.005 


0.100 


0.155 


W3-ALL 


0.010 


0.113 


-0.596 


-0.491 


-0.218 


-0.113 


0.008 


0.112 


0.181 


W3-RELAXED 


0.010 


0.110 


-0.731 


-0.525 


-0.203 


-0.109 


0.007 


0.108 


0.173 


W5-ALL 


0.009 


0.109 


-0.477 


-0.335 


-0.218 


-0.109 


0.004 


0.106 


0.185 


W5-RELAXED 


0.006 


0.103 


-0.451 


-0.287 


-0.198 


-0.101 


0.005 


0.100 


0.178 



